GRAVITY AND GRANULAR MATERIALS 

R P Behringer 1 , Daniel Hovell 1 , Lou Kondic 1 , Sarath Tennakoon 1 , Christian Veje 2 
'Department of Physics and Center for Nonlinear and Complex Systems, Duke University, Durham NC, 27708-0305 

bob@ phy.duke.edu 

^Center for Chaos and Turbulence Studies, Niels Bhor Institute, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 


We describe experiments that probe a number of 
different types of granular flow where either gravity 
is effectively eliminated or it is modulated in time. 
These experiments include the shaking of granular 
materials both vertically and horizontally, and the 
shearing of a 2D granular material. For the shaken 
system, we identify interesting dynamical phenom- 
ena and relate them to standard simple friction 
models. An interesting application of this set of 
experiments is to the mixing of dissimilar materials. 
For the sheared system we identify a new kind of 
dynamical phase transition. 


I. INTRODUCTION 

Granular Hows exhibit a rich phenomenology [1] 
that is still only partially understood. On the Earth, 
granular materials show dynamics that is often domi- 
nated by the effects of gravity. For instance, they are 
typically random densely packed at rest; if they initially 
have a large kinetic energy, so that they arc in a gas-like 
state, inelastic effects quickly dissipate that energy, 
leading back to the dense state. In a reduced gravity 
environment, the dynamics of granular materials is 
almost certainly different. The dense packed state is 
much less likely to occur, and other interesting phe- 
nomena are likely to be unmasked. For instance, in the 
present work, we have identified a second-order- like 
phase transition that is always masked by gravity for 
ordinary materials; this transition would be accessible 
in low gravity. Goldhirsch and Zanetti [2] have pre- 
dicted that as a granular gas cools, it will undergo 
another dynamical transition that will lead to spatial 
inhomogeneity, although gravity makes the study and 
interpretation of this effect on earth difficult. 

Understanding these phenomena would be scientifi- 
cally interesting, and quite possibly technically impor- 
tant. Such understanding would be crucial in any effort 
to mine in a reduced gravity environment, such as the 
Moon. 

In the experiments described below we have probed 
the effects of gravity in two ways. In the First, we have 
modulated gravity through shaking. In the second kind 
of experiment, we studied 2D granular systems where 
the “grains” slid on a smooth slippery surface that is 
horizontal. In this case, the surface carried the weight 
of the grains, and gravity played no role in setting the 
density of the grains. In the following two sections we 


highlight some of the recent results for each of these 
two types of experiments. 

II. MODULATION OF GRAVITY BY SHAKING 

Shaking provides a simple way of modulating grav- 
ity. Here, we vary gravity by shaking both in the verti- 
cal direction with a displacement 

z = A v cos(GV)> 0) 

and in the horizontal direction with a displacement 

x = Afccos(ov + <t>) • (2) 

The most important control parameters arc then the 
dimensionless accelerations Tj = A { CQ 3 /g, where i-h , 
o. 

For these experiments, we have developed a novel 
apparatus, Fig. 1, that allows us to provide independent 
2-axis shaking (a third axis could be added in future 
experiments). An additional feature of the apparatus is 
that it allows us to fluidize the grains by gas flow 
through a porous bottom plate. This fluidization is 
intended to reduce the contact forces between grains, 
and also reduces the effect of gravity. 



Linear beanngs 

Figure 1: Schematic of the apparatus. The rectangular 
cell is made of Plexiglas, and is mounted on a Plexiglas 
base of the same cross-sectional dimensions which is 
attached to a small table. The table is mounted on four 
linear bearings running on horizontal cylindrical guid- 
ance rods rigidly attached to a fixed bottom frame. The 
bottom section of the cell acts as a gas distributor. An 
electro-mechanical actuator, driven by a sinusoidal AC 
signal provides the horizontal driving. 
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Gravity plays a crucial role in granular friction, 
since gravity is responsible for the compaction of virtu- 
ally every earth-bound material. By modulating gravity, 
as in these experiments, we provide a useful probe of 
friction laws which date to the time of Coulomb, but are 
still only partially understood. The simple picture is 
that static friction provides a restraining shear force up 
to some maximum, Fj-~ after which slipping, i.e. 
failure can occur. At failure, the friction coefficient falls 
to its kinetic value, ji*. 

For purely horizontal shaking, we have studied [3] 
the onset of flow as T h T is increased from 0. Key find- 
ings of this work are contained in Figs. 2 to Fig. 5. The 
initial transition to flow is hysteretic, i.e. once flow has 
begun at T = T*, it can be sustained to lower T = T c . The 
hysteresis is lifted, however, if the layer is dilated about 
2% by gas fluidization. This amount of fluidization cor- 
responds to an effective decrease of the weight of the 
grains of about 40%. We interpret this result to mean 
that the primary source of hysteresis comes from the 
friction at grain contacts, with possibly a less important 
contribution due to the interlocking of grains. Here, a 
key point is that in ordinary descriptions of granular 
friction, it is not possible to distinguish between the 
normal forces generated by the interlocking of grains, 
and the tangenial friction forces at grain contacts. 
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Figure 2: Height, 177 of the fluidized layer vs. T^ for 
horizontal shaking. Data are for several materials, and 
show the hysteresis at onset. Flow begins at T* and 
ceases with decreasing T at T c . 
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Figure 3: Data for V* and T c vs. gas levitation pressure 
in units of the weight per area of the material. For 
modest amounts of gas levitation, the hysteresis in the 
transition to flow is removed. 


The overall flow pattern consists of sloshing motion 
which drives convection in the direction of shaking, 
plus a novel wall shearing effect that drives a weaker 
flow in the cross direction to shaking as well as down- 
ward at the sidewalls parallel to the shaking direction, 
Fig. 4. 



Figure 4: Sketch of time-averaged flow lines for hori- 
zontal shaking. The direction of shaking is the x-direc- 
tion. H is the height of the flowing (fluidized) layer. 
The strongest flow consists of sloshing of the grains in 
the x-direction, with down flow at the cndwalls, and 
corresponding inflow towards the center. There is also a 
shear-driven flow downward on average along the x-y 
bounding walls. 
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We have probed the latter effect through soft- 
particle MD simulations. These simulations are for 
7000 particles interacting via a soft-particle simulation 
[5]. The downward motion at the sidewall occurs 
because of dilation near that wall from shearing. Hence, 
near this wall, grains have an enhanced mobility, and 
can fall under gravity. We expect that this mechanism 
would occur routinely when there is a shear layer plus 
gravity. 



Figure 5: a) Experimental data (solid symbols) and 
MD simulations (open symbols) for shear-driven 
downward vertical velocity at the x-y walls vs. shaking 
velocity, A h (O h . b) Velocity profile, V y obtained from 
MD simulations, c) MD simulations show complex 
structure over a shaking period for the mean value of 
the particle momentum. 

For mixed vertical and horizontal shaking, a number 
of interesting effects occur [4], including the 
spontaneous formation of static heaps and at higher 
horizontal accelerations, the onset of sloshing flow. 
Static heaps form spontaneously at when < 

1, and T h is gradually increased from 0, Sloshing flow 
occurs at yet higher = r h2 when Coulomb friction 
can no longer sustain grains from sloshing up and down 
the inclined slope of the spontaneously formed heap. 

Below, we describe recent results in the context of 
Fig. 6. We understand the formation of static heaps at 
the point at which grains become unstable to the hori- 
zontal acceleration under the reduced effective gravity 
g(l - TJ that occurs as the shaker is accelerating 
downward. At the same time, if the shaker is driven 
horizontally in phase (0 = 0) with the vertical motion, 
then the grains will slide until they reach the appropri- 
ate end wall. This process stops when a steep enough 
incline has developed, so long as V h is not too large. 
However, if the V h is large enough, slipping will occur 
for either direction of horizontal motion, and sloshing 


of the grains results. Standard Coulomb friction pre- 
dicts all of these features qualitatively but not quantita- 
tively. Specifically, we consider a simple model 
consisting of a block resting on an plane inclined at an 
angle 0 and subject to uniform (i.e. non-oscillatory) 
accelerations of size gT v and gT^ in the vertical and 
horizontal directions respectively. The block will first 
slide on a horizontal surface when T/, = ^u(r c - T u ) (T c 
= 1), and we identify this event with the onset of static 
heap formation at T hl , Fig. 7. We identify sloshing in 
the granular system with a V h2 such that the block 
would slide for either direction of horizontal shaking, 
Fig. 8. A comparison of the model and the experiment, 
even when T c is allowed to be adjustable, shows that 
the model is only approximately correct. 



Figure 6: Data for the onset of flow under vertical and 
horizontal shaking. Here, is held constant at 0.68. 
Open symbols are pertain to T^ increasing from 0, 
starting from a flat heap. Solid symbols pertain to a 
similar increase in T^ starting from a heap inclined at 
the ordinary angle of repose. Solid curves are derived 
from simple Coulomb friction, based on the sketch in 
the inset. 

The shaker experiment provide several other 
interesting insights, although for space reasons, we do 
not provide extensive details here. One particularly 
interesting application of the shaking technique used 
here is for the mixing of dissimilar granular materials. 
Another system that we consider consists of a binary 
mixture of spheres that are identical except that some of 
the spheres have a higher coefficient of rolling friction 
than the other half. The spheres are then placed on a fiat 
surface that is shaken horizontally. We find interesting 
segregation effects that we will describe elsewhere. 
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Figurc 7: Comparison of data for the onset vs. F u 
for the spontaneous formation of a heap (T j) and the 
simple Coulomb friction model - i.e. a straight line. 
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Figure 8: Data for r«. the onset of sloshing, and the 
angle of the heap at T hl for = 0.68. Solid curves are 
based on the predictions of a simple friction model. 


The final experiment [6] consists of disks, a 2D 
granular material, that slide on a smooth slippery plane, 
and that are sheared in an annular geometry, Fig. 9. 
Specifically, the grains are sheared either by a rough 
inner wheel or by a surrounding rough outer ring. A key 
feature in these experiments is the absence of gravity, 
so that we can independently control the density. The 
disks are photoelastic, i.e. birefringent in proportion to 
the local deformation of the disks. We view them using 
a circular polariscope and use a novel technique to 


relate the photoelastic measurement to the applied 
force. In this technique, we exploit the fact that the 
transmitted light intensity, / depends locally on the 
local difference in the principle stresses with a disk. As 
the force at a contact increases, a series of light and 
dark bands appears within the disk, with the number of 
bands depending on the contact force. Hence, G, the 
gradient- square of 1YG = lAGi^, integrated over a grain 
size gives a good measure of the number of bands, and 
hence of the local applied force. In addition, each disk 
is marked with a dark bar so that we can track 
individual particle positions and orientations using 
video. From these measurements we deduce particle 
velocities, V and rotation rates (spins), 5. Here, we will 
focus on the velocity in the azimuthal direction, V 0 , 
where V 0 is the angular velocity of the grains. In 
addition, we measure the local packing fraction, y, 
which is effectively the density of disks. Here, we will 
focus on shearing by the inner wheel only, although it is 
also possible to shear by the outer ring. 


Figure 9: Sketch of 2D shear apparatus. Disks confined 
to a horizontal plane, and the region between wheel, A, 
and outer ring, B arc steadily sheared by the wheel. The 
wheel is roughened by teeth, E, and have small dark 
bars for tracking purposes. 


The figures below present some of the key results. 
Shearing by the inner wheel creates a shear band where 
the density (packing fraction, y) is significantly 
reduced, Fig. 10, and where there is a nearly 
exponential variation of the mean azimuthal velocity, 
Vq with distance, r/J, from the wheel, where d is the 
particle diameter (Fig. 11). In addition, the mean 
particle rotation rates, or spins, S show novel 
oscillations with distance from the wheel. There is a 
decreasing amount of dilation as y increases. We 
observe approximate but not perfect rate invariance in 
the V 0 and S, as we change the wheel rotation rate, Q, 
and we attribute the modest departures from rate 
invariance to small long time restructuring of the grains 
well away from the wheel. 
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Figure 10: Packing fraction profile, y vs. radial distance 
r/d ( d = particle diameter) from the shearing wheel, 
showing the formation of a shear band. 
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Figure 1 1 : Mean azimuthal velocity, V 0 , and mean spin, 
S vs. distance r from the shearing wheel in units of a 
disk diameter, D. Data are for different wheel rotation 
rates, Q = 2 x/x, where t is the time in seconds for one 
wheel revolution. Data are normalized to test for rate 
invariance, i.e., kinematic quantities such as S or V e 
should be linearly proportion to Q. 


Distributions for V 0 and S are also quite interesting, 
Fig. 12. These show a combination of stick and slip 
motion for the particles nearest the wheel. As the 
packing fraction y increase, there is transition from 
essentially complete slipping to nonslip with both 
rotation and translation of the disks. 
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Figure 12: Examples of the distributions for azimuthal 
velocity and spin of particles next to the shearing 
wheel. The packing fraction is y = 0.788 which is above 
the critical packing fraction. 


A particularly interesting feature of these 
experiments is the observation of a transition at a 
critical packing fraction, y c , that has some of the 
hallmarks of a second order phase transition. Here, y - 
y c is the order parameter. At y c r 0.77, the grains can 
just sustain shear, and the system is highly 
compressible. Also, there is critical slowing down: for 
instance the mean speed of the grains vanishes, Fig. 13. 
The distribution of forces changes qualitatively near y c . 
Well above y c , the distribution is qualitatively like the 
predictions of the q-model [7] of Coppersmith et al.: 
roughly exponential for large forces. Time scries for the 
stress and the corresponding spectra are similar to these 
that we found for sheared 3D materials [8]. 
Specifically, for a give shear rate and y, the spectrum 
typically varies as f' 2 at high frequencies, /, and as a 
powerlaw, P a/ “with 0 < a < 1 . Nearer the transition, 
the distribution is qualitatively different, and ongoing 
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measurements will be characterize the distribution in 
that case. We observe an interesting and apparently 
smooth change from slip to nonslip behavior as we pass 
beyond y c . 

Finally, if we define a slip event to be continuous 
decrease in the measured stress. Distributions of the 
size, M of slip events are exponentially distributed well 
above y c but appear to vary as AT 1 for small M when y 
is close to y c . 
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Figure 13: Mean azimuthal velocity and spin profiles 
for various packings fraction y showing slowing down 
near as y - y c from above. We estimate that y c = 0.77 
for this 2D system. 
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